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Turbulent curtains of smoke rise initially as flat plumes and, above a certain height, they become round 
plumes. The same evolution of cross-sectional shape is exhibited by jets issuing from flat nozzles. Here we 
predict based on principle that all such flows should evolve their cross- sectional shapes from flat to round 
(and not the other way) at a critical distance downstream, which is predictable. The principle is that the 
prevailing flow architecture provides greater access to the flow of momentum from the moving core (plume, 
jet) to the still surroundings. For turbulent plumes and jets, the transition distance scales with the long 
dimensions (L) of the two-dimensional (flat) heat sources and nozzles that drive them. For laminar jets, the 
transition distance scales with L Re, where Re is the Reynolds number based on nozzle velocity and the 
smaller dimension of the nozzle cross section. These predictions are confirmed by full numerical 
experiments of the three-dimensional flow fields of turbulent and laminar jets covering the Re range 10-10^. 

Organization in nature is in the eye of the observer. Where many see randomness and nondeterminism, 
others see organization and principle\ Nowhere is this more evident than in the discussion of tur- 
bulence. In spite of their obvious irregular appearance, turbulent flows are united by "organization" that 
makes their behavior expected, predictable, and reliable. 

For example, measurements of flow friction inside pipes and on solid surfaces have been lining up on curves 
that science and technology has used with confidence for more than a century. The large scale structure of jets and 
plumes, have become evident when the methods of measuring and visualizing turbulence improved. This 
empirical body serves as skeleton for the computational models and software packages used today to simulate 
turbulent flows at all scales, all the way to atmospheric circulation and the weather. 

The natural occurrence of organization is an integral part of the phenomenon of turbulence, and it is calling for 
a summarizing principle. Progress was made in this direction based on the constructal law of design and evolution 
in nature^"^, which holds that flow systems in nature, animate and inanimate, exhibit the natural tendency to 
evolve into successive configurations that offer greater flow access over time. The constructal law was used to 
predict and unify the turbulence phenomena, for example, the transition to turbulence, the size of the smallest 
eddy, the large scale structure of jets and plumes, and the correlations for friction and heat transfer. Reviews of 
these developments are available in Refs. 3-5. 

Here we draw attention to a pattern of turbulent flow evolution that went unquestioned, and shows how to 
predict it. The smoke from a row of closely spaced factory smoke stacks (or a brush fire) rises as a curtain, i.e. as a 
two-dimensional plume. It is only above a certain height that the curtain of smoke organizes itself into a round 
plume that, from a distance, looks like all the other plumes (e.g. Fig. 1). This observation is extremely common. 
Why is this? 

The reason is the natural tendency of flow systems to morph into configurations that facilitate the access to 
what flows. In jets and plumes, it is momentum that flows from the mover (the turbulent column) to the non- 
mover (the still surroundings). Momentum flows perpendicularly to the flow of fluid. This flow is called "mixing", 
or momentum transfer. When the lateral flow of momentum has greater access to the stationary environment, the 
fluid of the column mixes faster with the surrounding fluid, and its longitudinal speed decreases faster. The 
tendency of the flow is to morph its cross section such that the mixing becomes more effective, and the longit- 
udinal speed decreases faster. 

Turbulent free jets issuing from noncircular slots or nozzles have been studied because they are used in 
technical areas ranging from aerospace to chemical to mechanical engineering, such as aircraft design, gas turbine 
combustors, and cooling systems for turbine blades. For example, studies of free jets issuing from circular and 
noncircular nozzles documented the effect of the shape and size of the nozzle on jet development^"^. Other studies 
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Figure 1 | Above a certain height, all turbulent plumes have round cross 
sections. Top left: flat plume rising fi-om a row of smoke stacks. Top 
middle: round plume rising from a concentrated fire. Top right: plume 
above a brush fire. Bottom: The coalescence of plumes above the steam 
power plant at Duke University. (Photographs made by Adrian Bejan). 



documented the velocity and temperature fields by using linearized 
hot-wire anemometers^"^^. It is now estabhshed that there are three 
regions in the downstream direction of a three dimensional jet. The 
first is the potential core region, which contains the flow close to the 
jet exit with uniform longitudinal velocity, because mixing at the jet 
edges has not yet spread across the jet. Next is a region with decreas- 
ing longitudinal velocity that is dependent on nozzle configuration. 
This region bears the influence of the initial geometry of nozzle. 
Further downstream is an axisymmetric (round cross section) decay 
region, which contains axisymmetric flow that is independent of 
nozzle configuration^^"^^. Note the phenomenon that we identify, 
question and explain in this paper: the jet cross section evolves from 
nozzle-shaped to round, not the other way. 

At bottom, the phenomenon explained in this paper adds to the 
growing evidence that the tendency toward evolutionary organiza- 
tion in nature is governed by a physics principle. In particular, the 
tendency toward round cross sections in simple (one stream) con- 
figurations unites the animate with the inanimate fiow systems. 
Without exception, the cross sections tend to the round shape: blood 
vessels, pulmonary airways, subterranean rivers, large pores 
("pipes") in the hill slope, underground gafleries of moles and earth- 
worms, tunnels, intestines and many other features of animal design. 
In river channels, this natural tendency is made visible by the pro- 
portionality between depth and width in rivers of all sizes^'^^. 

Results 

Turbulent plumes. Plumes are vertical jets driven by buoyancy due 
to heating from below (Fig. 1) or the density difference between a 
lighter fiuid injected in a denser fiuid. The methods and results on 
plumes and jets are classical, well known, and not questioned in this 
paper. This is why they are relegated to the Supplement (for reviews 
see Refs. 5,20). Here the classical results are used in order to prove 
that the predicted evolutionary design from noncircular to circular 
fiow cross sections is correct. 



When the plume strength is large enough, the fiow is turbulent. If 
the turbulent field is time averaged^ the fiow occupies a mixing 
region the growth angle of which is approximately 20°. This rate of 
growth is a common feature of all turbulent mixing regions (shear 
layers, jets, plumes) 

A two-dimensional turbulent plume fills a wedge. Fig. 2. The key 
scale is the vertical speed of the plume (cf. Supplement), V2 ~ [q'gp/ 
(pcp)] where q' [W/m] is the strength of the heat source that drives 
the plume, g is the gravitational acceleration, and p, p and Cp are the 
coefficient of thermal expansion, density, and specific heat at con- 
stant pressure. 

A round turbulent plume (Fig. 2) fills a three-dimensional volume 
shaped as a cone rising above a point heat source (e.g. concentrated 
fire) of strength q[ W] . The vertical velocity (time averaged) along the 
plume centerline decreases with altitude (y) as V3 ~ [qgP/(pCpy)]^^^. 

The comparison between V2 and V3 is shown in Fig. 2. The fiow 
configuration that happens naturally is the configuration that offers 
better mixing with the ambient, and this means that the two-dimen- 
sional plume must evolve upward into a round plume. Below a crit- 
ical height, more effective mixing is offered by the two-dimensional 
plume, V2 < V3. Above that height, the more effective mixer is the 
round plume (V3 < V2). The height (ytr ~ L) where the transition 
occurs is predicted by intersecting V2 with V3. The transition height is 
of the same order as the larger length scale of the heat source that 
generates the two-dimensional plume, and does not depend on the 
plume strength q. 





Ytr - 



0 




Figure 2 | Turbulent plumes, flat (V2) vs round (V3). The competition for 
more effective mixing (the smallest v) dictates the transition from flat 
plume to round plume (and not from round to flat). 
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Turbulent Jets. Consider a two-dimensional jet issuing at speed Uq 
from a slit of spacing D2 and width L (Fig. 3). The scale of its 
longitudinal velocity decreases as x"^^^ in the longitudinal 
direction, cf Supplement, U2 ~ Uq [3y2D2/(4x)]^^^, with 72 = 7.67. 
If the nozzle is round with diameter D3 and velocity Uq, the 
longitudinal velocity decreases as namely U3 ~ Uq [y^^JsDs/ 
(4x)], with 73 = 15.2. 

The competition between U2 and U3 is shown in Fig. 3. Faster 
mixing occurs as the two-dimensional jet evolves into a round jet, 
not the other way around. The transition length [xtr ^ 7 3^3/ 
(47202)] is obtained by intersecting U2 with U3. The relation between 
D3 and D2 follows from the condition that the two jets must have the 
same strength at the nozzle, Uq (n/4) D3 = Uq D2 L, therefore 
Xtr ^ ylL/{ny2) = 9.6 L. The distance Xtr is proportional to the 
nozzle width L, and does not depend on the jet strength (D2Uq) or 
the nozzle Reynolds number. Re = U0D2/V. This conclusion is sim- 
ilar to the result for plumes. 

Laminar jets. The tendency to evolve from flat cross sections to 
round cross sections is also a characteristic of laminar jets, for 
which exact solutions exist in analytical form (cf. Supplement^^'^^), 
which are valid for Re > 1 when the jet region is slender and 
amenable to a boundary layer analysis. In the two-dimensional 
laminar jet the centerUne velocity decreases as x"^^^ in the 

downstream direction, U2 = 0.454 [UqD2L^/(vx)] where D2 
and L are the dimensions of the slit. Fig. 3. In a round jet issuing 
from a nozzle of diameter D3, the centerline velocity decreases as 
U3 = 3UqD3/(32vx). The intersection of the curves for U2 and U3 is 
qualitatively the same as in the lower graph of Fig. 3, and yields the 
transition distance Xtrjam = 0.132 U0D3/V =0.17 L U0D2/V, where 




0 Xtr X 



Figure 3 | Turbulent jets, flat (U2) and round (U3), and the transition from 
flat to round cross sections, in the direction of facilitating the momentum 
transfer from jet to ambient. 



we assumed that the two-dimensional jet has the same strength as the 
round jet. Unlike the transition lengths for turbulent flows, in the 
laminar jet the transition length Xtr/L increases with the nozzle 
Reynolds number. 

Discussion 

In summary, the empirical formulas available for plumes and jets 
were used to test and confirm the theoretical prediction that the flow 
configuration should evolve from the flat cross section to the round 
cross section, and not in the reverse direction over time. The trans- 
ition from one configuration to the next was determined by inter- 
secting the formulas available for the two types of cross-sectional 
shape, flat and round. 

The same predictions were further tested against numerical 
experiments with flat jets in the Re range 10-10^. With reference 
to the configuration shown at the top of Fig. 3, the jets issued from 
a slit with the ratio L/D2 = 10. The numerical formulation and 
software are described in Supplement. The simulations produced 
time-averaged velocity distributions through the entire flow field. 

For each Re, the objective of one simulation was to document the 
evolution of the shape of the jet cross section with longitudinal posi- 
tion (x, collinear with U2 in Fig. 3) downstream from the nozzle. The 
cross section of the jet was defined in a constant-x cut as the area (A) 
enclosed by the loop along which the longitudinal velocity was equal 
to 1/2 of the centerline velocity. Samples of this conventional cross 
section are shown for Re = 500 in the upper frame of Fig. 4. Later, 
this conventional velocity was set equal to 1/4 of the centerline velo- 
city, in order to verify that this convention has no effect on the 
conclusions regarding the transition from flat jet to round jet. 

The evolution of the shape of A was monitored by plotting the 
shape parameter defined here as p/A^^^, where p is the perimeter of A. 
When the cross section is flat, p/A^^^ is large, for example p/A^^^ = 22/ 
-j^Qi/2 ^ y ^Yie plane of the nozzle. When the cross sections are 
round, the shape parameter is small, p/A^'^ = 2n^'^ = 3.5. 

Figure 4 confirms the two parts of the theory. The transition from 
flat cross sections to shapes that are not flat occurs when x/L is of 
order 10. This is in accord with the theory, Xtr ~ 9.6L. The shape in 
the transition range is not perfectly round. It oscillates once or twice 
between a horizontal ellipse and a vertical ellipse, however, as x/L 
increases the elliptical shape becomes progressively rounder. The 
second theoretical aspect is that the Xt^/L value should not depend 
on Re. This is also confirmed by the variation of the shape parameter 
p/A^^^ versus x/L in the Re range 500-10^. 

Further support for the theory is provided by the lower graph of 
Fig. 4, which shows the corresponding results for laminar jets in the 
range Re = 10-30. The transition from flat to round jets is in accord 
with the theory, Xtrjam ~ 0.17 Re. For example, the predicted Xtr/L for 
Re = 10 is of order 1, in accord with the x/L location of the elbow in 
the Re = 10 curve. For laminar jets, the predicted x^^/L should 
increase in proportion with Re. This is also confirmed by the results 
of Fig. 4. The elbows of the curves migrate to the right by a factor of 3 
on the logarithmic abscissa, while the Re values increase by the same 
factor. 

Conclusion 

In this paper we used a common example of natural organization in 
order to demystify the phenomenon of evolution in nature. The 
example was inanimate: the evolution of all plumes and jets in one 
direction in time, from any cross-sectional shapes to round cross 
sections. The arrow of time points in one direction, and this direction 
is captured by the constructal law of design evolution in nature. 

The choice to illustrate evolution with an inanimate flow was 
intentional, because the prevalent view is that evolution is a phenom- 
enon of biology, and that it cannot be witnessed in our lifetime. We 
showed that in plumes and jets we see evolution in real time. We also 
see the principle that allows us to predict turbulence as a phenom- 
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Figure 4 | The evolution of the cross-sectional shapes of turbulent jets in the range Re = 500-10^ and laminar jets in the range Re = 10-30. 



enon of physics. The constructal law was used earUer^ to predict the 
laminar-turbulent flow transition, the transition from splat to splash 
as a liquid droplet hits a wall, the dendritic architecture of the snow- 
flake, river basins and lightning trees. It can be used to predict the 
structure of two -phase jets^^'^^. 

In sum, evolution is a universal phenomenon that unites the anim- 
ate and inanimate flow designs of nature. Evolution unites biology 
with physics. In the human realm of life (movement) on the globe, 
this universal phenomenon is visible as technology evolution, science 
evolution, language, writing, the rule of law (government), and civil- 
ization. All forms of evolving organization obey the arrow of time^^ 
illustrated in this paper. 
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